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ABSTRACT 

A comparison between linear stability analysis and observations of pulsation modes in five 
5 Scuti stars, belon ging to the same cluster, is presented. The study is based on the work by 
iMichel et all (Il999l) . in which such a comparison was performed for a representative set 
of model solutions obtained independently for each individual star considered. In this 
paper we revisit the work by Michel et al. ( 1999) following, however, a new approach 
which consists in the search for a single, complete, and coherent solution for all the 
selected stars, in order to constrain and test the assumed physics describing these ob- 
jects. To do so, refined descriptions for the elfects of rotation on the determination of the 
global stellar parameters and on the adiabatic oscillation frequency computations are 
used. In addition, a crude attempt is made to study the role of rotation on the prediction 
of mode instabilities. T he present results are found to be comparable with those reported 
bv lMichel et a l.' ('1999'). Within the temperature range log Teff = 3.87-3.88 agreement be- 
tween observations and model computations of unstable modes is restricted to values for 
the mixing-length parameter (Tnl - 1 -50. This indicates that for these stars a smaller 
value for o-nl is required than suggested from a caUbrated solar model. We stress the 
point that the linear stability analysis used in this work still assumes stellar models without 
rotation and that further developments are required for a proper description of the interaction 
between rotation and pulsation dynamics. 
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1 INTRODUCTION 

With spectral types from A2 to FO, the 6 Scuti stars are pulsating 
variables both in the main sequence and in the sub-giant evolution 
phase (hydrogen shell burning). They are considered as good can- 
didates for asteroseismic studies. Such stars show modes excited 
over a wide frequency range, including mixed modes known to be 
sensitive to the structure of the deep interior. These characteristics, 
together with the absence of magnetic field or metallicity peculiari- 
ties make them suitable for the asteroseismic study of hydrodynam- 
ical processes occurring in stellar interiors, such as convective over- 
shooting, m ixing of che mical species, and redistribution of angular 
momentum I Zahnll992h . Recently, Suarez, Goupil & Morel (2006) 
reported that the efi^ect of the distribution of angular momentum 
on adiabatic oscillation frequencies in 6 Scuti stars is significant 



and that it can be detected by space missions such as CoRoT. Due 
to the complexity of the oscillation spectra in 6 Scuti stars, their 
pulsating behaviour is not completely understood (see 'Cox"200^ 
for a complete review of unsolved problems in stellar pulsation 
physics). In particular, great efforts are made nowadays to solve the 
problem of mode identification. Ground based multisite campaigns 
are regularly o rganised withi n coordinated networks, e.g.: STEPHI 
(Michel et al.loOO) or DSN (Bre2er"2000';'Handler"2000^. Only a 
few tens of the low-degree modes have been detected from ground- 
based observation s, with a maximum of around 75 modes for the S 
Scuti star FG Vir tereger et al.'2005'). The recently launched space 
mission (Baglin et al. 2002) represents a unique opportu- 

nity for investigating such stars considering its much lower detec- 
tion threshold and for obtaining data from quasi-uninterrupted time 
series over 5 months. 
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Table 1. Detected frequency peaks in Praesepe target stars. 
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Additional uncertainties arise from the effect of rapid rotation, 
both directly on the hydrostatic balance in the star and through mix- 
ing caused by circulation or instabilities induced by rotation. The S 
Scuti stars are commonly fast rotators (100 < Vsini < 200 km s"'), 
and consequently the symmetry of the multiplets is broken by the 
rotation. In the framework of a perturbation analysis, the second 
order effects induce strong asymmetries in the splitting of mul- 
tiplets Isai o 1981; Dziembowski & Goode 1992) an d frequency 
shifts must not be neg lected even for ra dial modes jSoufi et all 
Il995h . In this context, iMichel et al. I ( fT999i) , hereafter M99, pro- 
posed a technique to estimate and correct the effect of fast rota- 
tion on the determination of fundamental parameters for stars in 
clusters. Applying this technique to a set of 6 Scuti stars belonging 
to the Praesepe cluster, the authors showed that it was possible to 
reach a reasonable agreement between ranges of observed oscilla- 
tion modes and overstable radial modes predicted by a linear stabil- 
ity analysis. However, agreement was found only for certain values 
of ffNL, the mixing-length parameter that was used in the non-local 
time-dependent convection treatment of the stability computations. 
This implies that models for convective heat transport in the outer 
stellar layers may be calibrated against observations. 

In M99, different series of stellar model solutions were ob- 
tained, which assumed a wide range of physical parameters, 
e.g. overshooting, mixing-length parameter, metalHcity, stellar 
age, etc. Considering these results as our reference domain of 
possible solutions, the present work revisits the work by M99 
but adopts an improved approach. The main idea is to search for 
one particular solution that explains the whole set of observations, 
instead of sets of individual solutions for each star. To do so, refined 
techniques for modelling intermediate mass stars are required, in 
particular such techniques that take different effects of rotation into 
account. In concrete terms, in M99 the photometric parame- 
te rs are corrected for the effect of rotation using the results 
of iMaeder & PevtremannI jl995l) who applied the von Zeipel 
gravity-darkening law to estimate the emergent flux of a ro- 



tat ing star. In the present work w e follow the method described 
by IPerez Hernandez et all jl999l) (hereafter PII99) which im- 
proves such a calculation by means of, among other modelling 
aspects, up dated m odels and consider the gravity-darkening 
law given bv lClaret A . ( 1998), which, unlike the von Zeipel law, 
is also valid for non-fully radiative stars. Moreover, in M99 non- 
rotating equilibrium models were used whereas here we employ 
pseudo-rotating models. Such pseudo-models take first-order 
effects of rotation into account by means of an effective gravity 
(see Sect. l3.2t . One of the major improvements with respect to 
M99's work lies in the computation of adiabatic oscillation fre- 
quencies. In M99 the effects of rotation on the oscillation fre- 
quencies (up to the second order) are included in the manner 
o f Perez Hernandez et al»tI995.) . considering the perturbed fre- 
quency v'„ as 

y'„ = y„ + (A„ + B)^ , (1) 
v„ 

where v„ is the unperturbed frequency, Vpo, is the rotation 
frequency, i? is a constant with an asymptotic value of 5/3, 
and A„ corresponds to a coefficient depending on the eigen- 
functions, and thus, on the stellar model. These coeffic ients 
are o btained by interpolating in the values computed by ISaiol 
for an n = 3 polytrope. In the present work, following 
Ipziembowski & Goode ( 1992) and Soufi et al. ( 1998) (see also 
Suarez, Goupil & Morel 2006) a complete treatment of second 
order effects (including near degeneracy) is used. Furthermore, 
this formalism takes the effect of the star deformation due to 
rotation into account. 

With these improvements a similar methodology than that 
adopted in M99 is used here to study five 6 Scuti stars of the 
Praesepe cluster, four of which have already been included in 
the sample considered by M99. The fundamental parameters 
of these stars are determined from taking into account the ef- 
fect of rotation on the photometric observables using the up- 
dated method by PH99. These fundamental parameters are 
then used to build representative asteroseismic models for each 
star, consisting of pseudo-rotating equilibrium models (instead 
of non-rotating equilibrium models used in M99) and their cor- 
responding adiabatic oscillation spectra (properly corrected for 
the effect of rotation). These asteroseismic models are then used 
to determine the observed ranges radial orders, which, as it was 
done in M99, are then confronted with mode instability predic- 
tions obtained from a linear stability analysis. We note that, as 
in M99, the instability predictions are carried out using equiv- 
alent envelope models which do, however, not take the effect of 
rotation into account. This is so because, up to date, there are no 
reliable theories available which describe the effect of rotation 
on mode stability. Nevertheless, a crude estimate of this effect is 
addressed, which assumes that mode stability depends predom- 
inantly on the effective temperature of the model ( Pamiatnvk^ 
Il975[) . 

The paper is structured as follows: In Section[2]the main char- 
acteristics of the Praesepe cluster and the observational material 
are presented. Equilibrium models and isochrone computations are 
discussed in Section [3] Different aspects of the oscillation compu- 
tations are described in Section |4] In Section |5] the consistency of 
our solutions with the results obtained by M99 is examined. Then, 
a detailed discussion of the present results is given. Finally, conclu- 
sions and perspectives are summarised in Section|6] 
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Table 2. Photometric data for the five selected S Scuti stars. The columns 
represent name, HD number, absolute magnitude (in mag), Bj - V\ colour 
index (in mag), and the projected velocity Vsin; (in kms"') of the star. 



Star 


HD 




B2'Vi 


Vsinr 


BWCnc 


HD 73576 


2.18 


0.063 


170 


BSCnc 


HD 74763 


2.19 


0.048 


135 


BVCnc 


HD 73798 


2.36 


0.079 


110 


BUCnc 


HD 73450 


1.37 


-0.004 


205 


BNCnc 


HD 73746 


1.52 


0.016 


200 



3 MODELLING 

In the process of searching representative models for the five se- 
lected 5 Scuti stars, we take advantage of their cluster membership 
in several aspects: firstly, the chemical composition and age are as- 
sumed to be the same (to a relatively good extent) for all objects. In 
addition, metallicities and distances are commonly estimated better 
for stars in clusters than for field stars. Secondly, the photometric 
parameters of target stars are corrected for the effect of rotation 
following M99 and PH99. Such corrections also provide improved 
estimates for the stellar parameters required for the modelling. 



3.1 Stellar models and isochrones 



Praesepe 




-0.15 -0.10 -0.05 -0.00 0.05 0.10 0.15 
B2-V^(inag) 

Figure 1. Colour - magnitude diagram for the Praesepe cluster. Crosses rep- 
resent the observations of stars with confirmed membership. Filled circles 
correspond to the 6 Scuti stars selected for the present work. All objects are 
corrected for distance. Curves represent isochrones of relevant ages com- 
puted from models with ffMLT = 1.614 and rfoy = 0.2. 



2 OBSERVATIONAL DATA 

The selection of stars considered here contains five 6 Scuti stars 
belonging to the Praesepe cluster. Four of these stars: BWCnc, 
BSCnc, BUCnc and BNCnc (HD 73798, HD 73450, HD 73576, 
and HD 73763 respecti vely) were observe d by several campaigns 
of the STEPHI network jMichel et al.ll995il The fif th star, B V Cnc 
(HD 73746), was observed by Frandsen et alj ( l200lh . The observed 
frequencies of these stars are listed in Table[T] 

Regarding the cluster parameters, its dis tance modulus (6.28 + 
0.12 mag) was taken from iRobichon et al. 1 l i 19991) who derived it 
from mean cluster parallaxes computed from Hipparcos intermedi- 
ate data. We assume a metallicity value [M///] of 0.170 which was 
derived from Geneva photometry of single stars with spectral types 
in the range F4-K3 (Grenon 2000; Robichon et al. 1999). Apparent 
magnitudes and indices were taken from the Rufener's catalogue 
( lRufenerlll988l) . In Table [2] the specific photometric data for the 
five selected 6 Scuti sta rs are listed. The p rojected velocity values 
(Vsini) were taken from lRover et alj ( |2002|) . 



To characterise theoretically the observed stars within the clus- 
ter we first construct standard non-rotating models and isochrones 
that take into account the observed parameters presented in Sec- 
tion (2] and second, we also build models including uniform ro- 
tation (pseudo-rotating models) as described in Section [T2l The 
evolutionary stellar models are computed with the CESAM code 
jMorell 1997). Opacity table s are taken from the OPAL pack- 
age M 

esias & Roger sl l 19961) , complement ed at low temperatures 
(TsJ lO'' K) by the tables provided by lAlexander & FergusonI 
[1994'). The atmosphere is constructed from a grey Eddington T - t 
relation. The Praesepe's metallicity value of ZjX is derived from 
the [M///] value given i n Section |2l assuming (Z/X)o = 0.0245 
( iGrevesse & Noelsil 19931) . Y^, = 0.235, and Zp, = for the he- 
lium and heavy element primordial mass abundances, and the value 
AK/AZ = 2 for the enrichment ratio. This yields an initial helium 
abundance Y = 0.285 and a heavy element abundance Z = 0.025 
for the cluster. Convection is treated with a local mixing-length 
model teohm- Vitense" 1958^ and we assume various values for the 
mixing-length parameter q-mlt = II H^, where / is the mixing length 
and Hp is the pressure scale height. In particular, values between 
0.5 and 2, including the calibrated solar value of aMLx = 1.614, 
are considered. For the overshoot parameter rfov = low I (low be- 
ing the penetration length of the convective elements) we assume 
values between 0.0 and 0.3. 

For the construction of non-rotating isochrones, various sets of 
evolutionary sequences were computed for masses between 1 Mg 
and 5 Mq, evolved from the zero-age main sequence (ZAMS) to the 
sub-giant branch. Isochrones were then obtained using the Geneva 
Isochrone Code. Three representative isochrones around the age of 
the Praesepe cluster are depicted in Fig.[T] after having previously 
transformed their effective temperatures and luminosities into the 
Geneva photometric system. Th is transformation is p erformed us- 
ing the cali bration pro posed by Schmidt- Kaleil i ll982h for My, and 
by KunzherZI ( ll997b for (S2 - Vi). Potential binarity and the ef- 
fect of fast rotation, discussed in Section [T2l are expected to in- 
duce systematic shifts of the models towards higher luminosities 
and lower effective temperatures when compared with single non- 
rotating stars. Therefore, in order to avoid both effects, the fit of 
isochrones has to be made by adjusting the isochrones to the bottom 
envelope of the cluster in the colour-magnitude diagrams (hereafter 
CM diagrams). 

Considering the physical parametrisation described above, a 
range of ages between 650 Myr and 750 Myr is then found to be 
representative of the cluster (Fig. [T). This range is in agreement 
with values found earlier in the literature for the Praesepe cluster. 
The age will be adjusted more carefully in Section [T2l 
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Figure 2. Colour - magnitude diagram sliowing the corrections (two seg- 
ments per star) of tlie observed photometric parameters for the effect of ro- 
tation. Circles represent the observed 6 Scuti stars. Dashed and continuous 
segments correspond to empty and filled circles, respectively. Tiie three 
pairs of segments on bottom are associated with BS Cnc, BW One, and 
BV Cnc. The two remaining ones on top are associated with BU Cnc 
and BN Cnc. Thicli segments correspond to corrections obtained from 
an increase of Vsin! by -1-10%, and the thin segments correspond to cor- 
rections obtained from decreasing Vsin; by -10%. The continuous curve 
corresponds to a 650 Myr isochrone (see Fig.[T). 



3.2 Stellar parameters corrected for rotation effects 

The modelling of the observed 6 Scuti stars requires the knowl- 
edge of certain basic stellar parameters including the values for 
mass, effective temperature and luminosity (and/or gravity). Two 
of these quantities can be determined from high resolution spec- 
troscopy or, as in the present case, by means of photometric ob- 
servations. However, as shown in M99, fast rotation must be taken 
into account when estimating fundamental parameters from photo- 
metric indices. In order to take this into account, the authors also 
proposed a method for correcting the photometric indices for the 
effect of fast rotation. Here we use an improved version of this 
method, developed by PH99 who considered three main effects: 
Firstly, the shape distortion caused by rotation is such that a 
star seen pole-on has a larger projected area and hence looks 
brighter than if it were seen equator-on. Secondly, the flux (and 
hence the effective temperature) is larger at the poles than at 
the equator (gravity-darkening effect), and thirdly, rotation de- 
creases the intrinsic luminosity and increases the mean radius 
of the star. As a result of these effects the location of a rotating 
star in a colour-magnitude diagram (and thereby in an HR di- 
agram) depends not only on the angular velocity of the star but 
also on its angle of inclination ;'. The basic idea of the method is 
to obtain, for a given rotating star and for a given Vsini value, a 
segment covering the different potential positions in the colour- 
magnitude diagram (i.e. for the corresponding putative stars 
which are not rotating). The corrections so obtained for the se- 
lected 6 Scuti stars are displayed in Fig.|2l in which the positions 
of the non-rotating co-partners are varied with the angle of in- 
clination 7, and with the rotation rate ai, defined as oj = Q./Q.^, 
where Q is the angular rotational velocity of the star, and Sl^ 



is its break-up rotational velocity^ Although the corrections 
are model-dependent, for a given isochrone, this dependence 
is weak and the uncertainty introduced smaller than that com- 
ing from the observational error in Vsin/. An error on the Vsin/ 
values of ±10% is assumed, which yields two segments per star. 

After an iterative process (more details in PH99) it is possible, 
in most of the cases, to converge to an intersection between the 
isochrone and the segments. The non-rotating co-partners of our 
observed stars closest to the intersection between the isochrone and 
the segments give us an estimate of both / and a>. Additionally, 
the isochrone models at the intersection provide estimates of the 
required remaining stellar parameters. 

As mentioned in PH99, the photometric correction for rotation 
depends on the selected isochrone. In the present case, the convec- 
tion parameters ttMLT and d„y were adjusted within the ranges given 
in the previous section; the best solution (given by this correction) 
was obtained for qtmlt = 1.614 and d„y = 0.2, leading to an optimal 
age of 650 Myr for the cluster. This solution was selected under 
two main criteria: first, a standard fitting of isochrones, and sec- 
ond, the existence of suitable photometric corrections (see PH99 
and Suarez et al. 2002 for more details) for the whole sample of 
stars. Both criteria make the solution converge rapidly to an age of 
650 Myr (±20 - 40 Myr). The age uncertainty of 20-40 Myr can 
be neglected in terms of global characteristics of the non-rotating 
co-partners (see the influence of ag e on photometric co rrections for 
rotation in different open clusters in lSuarez et al 12002). For this so- 
lution, the following range of masses are found: from M = 2.05 - 
2.1 IMo (om high mass stars) to 1.64- 1.65M0 (our low mass stars), 
which is within the range of masses obtained in M99. In particular, 
the high mass models are found to be representative of BU Cnc 
& BNCnc; the M = 1.70- 1.73 Mq models are representative of 
BW Cnc & BS Cnc, and finally, the M = 1 .64 - 1 .65 Mq models are 
representative of BVCnc. The characteristics of the different mod- 
els are summarised in Table [3] We recall that the different ranges 
listed in Table [3] correspond to solutions in which Vsin; was var- 
ied by ± 10% from its observed value. As explained in PH99 and 
ISuarez et al. I l l2002h , this uncertainty largely dominates the errors 
of the method. 

We then compute stellar models including uniform rotation 
for the five considered stars, adopting the corrected values for their 
masses and ages obtained from the method described above. To 
take first order effects of rotation into account, the equilibrium 
equations are modified in the CESAM code in the manner de- 
scribed in iKippenhahn & Weigerll ll99&). In particular, the spher- 
ically averaged contribution of the centrifugal acceleration is in- 
cluded by means of an effective gravity geff = g - ■^c(r), where 
g is the local gravity, r is the radius, and ^^(r) = | rQ}(r) is the 
centrifugal acceleration of matter elements. This spherically aver- 
aged component of the centrifugal acceleration does not change the 
order of the hydrostatic equilibrium equations. From now on we 
shall refer to such models as 'pseudo-rotating' models. Although 
the non-spherical components of the centrifugal acceleration are 
not considered in the equilibrium models, they are included in the 
adiabatic oscillation computations by means of a linear perturba- 
tion analysis according to Soufi, Goupil & Dziembowski (1998) 
(see also Suarez, Goupil & Morel 2006). 



^ This corresponds to the rotational velocity that the star would have if the 
centrifugal force is balanced by the gravitational attraction at the equator. 
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Table 3. Characteristics of the non-rotating counterparts for the selected 6 Scuti stars of the Praesepe cluster. These characteristics are obtained with the method 
for correcting the effect of rotation on photometric parameters (Section [3.2t . The indicated ranges correspond to solutions, from left to right, for which Vsin/ 
was varied by +10% and -10% from its observed value (see more details in Section[T2). The best results with this method are obtained for an age of 650 Myr 
for the cluster. The different columns are: the star identification from the General Catalogue of Variable Stars (GCVS), the absolute magnitude My; the colour 
index Bt — Vi ; the mass M (in units of Mq); the radius R (in units of Rq); the effective temperature Teff (log, in K); the gravity g (log, in cm s^); the rotation 
rate a» (in units of the break-up frequency n^); the cyclic rotation frequency Vpot (in juHz), and finally the inclination angle i of the star. 



ID Mv 52-^^1 M R Tefi g oj Vyo, i 



BUCnc 


1.29- 


-1.36 


-0.070- 


(-0.057) 


2.09 


-2.11 


2.52- 


2.56 


3.91 


-3.91 


3.95 


-3.94 


0.95- 


0.95 


19.05 


-14.70 


56.27 


-51.17 


BNCnc 


1.36- 


-1.35 


-0.062- 


(-0.055) 


2.05 


-2.05 


2.42- 


2.42 


3.91 


-3.91 


3.98 


-3.98 


0.90- 


0.80 


18.93 


-16.83 


90- 


66.30 


BWCnc 


2.19- 


-2.23 


0.011 


-0.014 


1.73 


-1.70 


1.87- 


1.84 


3.88 


-3.88 


4.13 


-4.14 


0.90- 


-0.90 


24.14 


-24.71 


42.63 


-48.25 


BSCnc 


2.21- 


-2.24 


0.011 


-0.014 


1.73 


-1.70 


1.86- 


1.83 


3.88 


-3.88 


4.13 


-4.14 


0.80- 


-0.80 


20.04 


-20.35 


45.47 


-40.22 


BVCnc 


2.41- 


-2.42 


0.043 


-0.044 


1.65 


-1.64 


1.75- 


1.74 


3.87 


-3.87 


4.16 


-4.17 


0.80- 


-0.80 


24.58 


-24.72 


34.70 


-31.42 
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Figure 3. Observed and predicted (using linear stability analysis) ranges of unstable radial modes for the selected 6 Scuti stars (« i is the lowest value, and nj 
the largest value of the radial order of the unstable modes). In panel (a) we re-plot the results obtained by M99 for BU Cnc, BN Cnc, BS Cnc and BW Cnc, and 
in panel (b) the current results are displayed. Filled circles represent the observed ranges. Rhombus and squares correspond to predicted radial order ranges 
foraNL = 1.89 and Q-NL = 1.50, respectively. Each diagonal-dashed Une represents the width (in radial orders) of the represented ranges. 



4 DETERMINATION OF OBSERVED AND 

THEORETICAL RANGES OF UNSTABLE MODES 

We compute synthetic adiabatic oscillation spectra correspond- 
ing to the pseudo-rotating evolutionary models described in Sec- 
tion |3.2| We restrict our computations to modes with degree £ ^ 2 
since modes with £ ^ 3 are generally considered to be invisi- 
ble in photometry (only visible in spectroscopy) for distant stars 
due to geometric cancellation. The ranges of observed radial and 
nonradial modes are similar to the predicted ranges of computed 



overstable radial modes. This is to be expected, because driving 
and damping in <5 Scuti stars takes place predominantly in the 
Hell ionisation zone, which is rather close to the stellar surface, 
where the vertical scale is much less than the horizontal scale 
of the oscillations and when £ is low the modal inertia is quite 
insensitive to degree {. It is therefore plausible to use only ra- 
dial modes in the computations the res ults of which are also ap- 
plicable to modes of low degrees ( eg. iDziembowski et al]|200ll : 
iDaszvnska-Daszkiewicz et"al]l2005[) . Then, proceeding as in M99, 
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for each star, the range of predicted unstable radial modes (accord- 
ing to linear growth rates) is compared with observed mode ranges. 

4.1 Determination of observed ranges of radial modes 

Theoretical adiabatic oscillati on spectra are computed with the os- 
cillation code FiLou ( Tran Minh & Leorj 11995 ; Suarez 2002). As 
mentioned in Section[l] we are using an improved version of M99's 
treatment for second-order effects of rotation on the adiabatic os- 
cillation frequencies. In particular, we use th e complete treat- 
ment of second-order effects of rotation by Suarez et al. ( 2006), 
which is based on the form alisms by Dziembowski & Goode 
jl992h and lSoufl et ai] jl998h . Furthermore, the theoretical adi- 
abatic oscillation frequencies are computed from the pseudo- 
rotating equilibrium models (Section |3.2> , whereas M99 con- 
sidered only non-rotating equilibrium models. We recall that 
these pseudo-rotating equilibrium models, representative of 
each star, are computed from the rotationally corrected stellar 
parameters of Sect. 13.21 

For each star the range of radial orders of the observed modes 
is then determined in the manner of M99, i.e. the locations of the 
peaks in the observed power spectrum are compared with the loca- 
tions of radial modes in the theoretical eigenspectrum of the corre- 
sponding pseudo-rotating models. 

4.2 Predicted unstable modes 

The stabil i ty co mputations are carried out in the manner of 
lBalmforthl ( ll992l) . We include the Lagrangian perturbations of the 
turbulent fluxes, i.e. the convective heat flux (5Fc and turbulent 
pressur^ (5p, in the stability computations acc ording to the n on- 
local, time-dependent mixing-length model by iGouglJ ( fT977al lbl). 
The equilibrium envelope models are obtained by specifying mass, 
luminosity and effective temperature as provided by the evolution- 
ary computations of Section 3, but they are not corrected for 
rotational effects. The radiation field is treated in the Eddington 
approximation, and the atmosphere is considered to be grey and 
plane-parallel. We adopt the nomenclature of M99, i.e. q-nl = IjH^ 
is the mixing-length parameter of the non-local convection model 
used in the stability computations. The non-local mixing-length 
parameter c^nl is calibrated to the same depth of the outer con- 
vection zone as suggested by the evolutionary computations (see 
Section [3) wh i ch use the standard mixing-length formulation by 
lB6hm-Viten"s3 ( ll958h and a local mixing-length parameter q-mlt = 
1.614. The so obtained calibrated value of q-nl = 1.89. To study the 
effect of varying the mixing length on mode stability we computed 
a second series of stellar models with ttNL = 1.50. Further detail s 
on the stability computations can be found in lHoudek et al. Il ll999l) . 



5 COMPARISON BETWEEN OBSERVED AND 
PREDICTED RANGES OF UNSTABLE RADIAL 
MODES 

We proceed as in M99 and compare ranges of observed and pre- 
dicted radial orders n of unstable modes for two values of a^i^: 1.89 
and 1.50 (see Section |42t . We recall that this comparison does not 

^ Pi = pw^ h the (r, r)-component of the Reynolds stress, where p is the 
density of the mean stratification and w is the rms vertical component of the 
turbulent velocity field. 



imply that all observed modes are identified as radial modes, but 
their frequency ranges are analysed in terms of unstable modes de- 
scribed within a range of radial modes (for details see M99). 

The errors in the calibration from colour indices to effective 
temperatures can reach 150 K. It is therefore adequate to accept an 
uncertainty of ± 1 w in the determination of the range of radial or- 
ders. In addition, this error also accounts for (1) the difference of 
ranges estimated in the co-rotating and the observer's frame, i.e., 
the rotational splitting; for modes up to f = 2, this represents a fre- 
quency range shift up to v^o, C ~ 40 pHz, with Vj^, being the stellar 
rotational frequency (up to ~ 20 pHz); and (2), for the possible 
C = 1 and 2 modes with frequencies close to the boundaries of 
the radial order ranges. 



5.1 Consistency with the M99 results 

As a first step, we examine the consistency of the present results 
with the results by M99. We do so by comparing the ranges of radial 
orders (observed and predicted ranges) obtained here (Fig.[3l5) with 
the ranges reported by M99 (Fig.[3}i) for BU Cnc, BN Cnc, BW Cnc 
and BS Cnc . In both panels, the observed ranges of radial orders 
are represented by filled circles; those obtained from the stability 
analysis using q-nl = 1.89 and 1.50 are illustrated by rhombus and 
squares, respectively. The coordinates (ni, «2) correspond to ranges 
of radial orders between iii and 112. 

When performing this comparison, it has to be kept in mind 
that the present work includes various improvements in the stellar 
modelling over M99 and consequently represents a more specific 
and more precise solution, whereas the M99 results represent an 
extended and conservative set of solutions in terms of model pa- 
rameters. In this view, we are interested to check if our new solu- 
tions do belong to the set of solutions by M99 (and if so, to what 
extent), or if they differ significantly. 

Considering the observed ranges. Fig. [3] shows that the new 
solutions lie within those reported by M99. The only significant 
difference is observed for BN Cnc and which has to be attributed to 
the updated (high) Vsini value adopted in the present work. 

As for the instability predictions, our results for the two con- 
sidered values of ttNL are in general similar, within ±ln, to the 
results reported by M99. 

5.2 Discussion of the present results 

Focusing on the present results, we come to the comparison be- 
tween observed and predicted ranges which is illustrated in Fig.[3]5. 
Generally, the observed ranges are in good agreement with the the- 
oretical predictions using q-nl = 1.50, whereas with q-nl = 1-89 
the predicted ranges are in grave disagreement with the observa- 
tions for the intermediate mass stars BW Cnc and BS Cnc. In M99 
this agreement (and disagreement) was observed for the set of four 
stars. However, in M99 the comparison between observations and 
theory was performed for a representative set of model solutions 
which were obtained independently for each individual star. On 
one hand, M99 reported that the observed and theoretical results 
could not be made to agree with q-nl = 1.89. On the other hand, 
their comparison was in reasonable agreement for models com- 
puted with Q-NL = 1.50 which, however, does not necessarily mean 
that they found one consistent model solution for all four objects 
(stars) all of which having the same age and metallicity. The results 
reported here, particularly the fact that we also do not find agree- 
ment between observations and theory for models obtained with 
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Figure 4. Same as Fig. [3] but for models with eflfective temperatures ob- 
tained from photospheric parameters that are corrected for rotational effects 
(more details in text). 



ttNL = 1.89, support the conclusions reported by M99, but lim- 
ited to only one consistent model solution for all considered stars, 
which is the main objective of the present work. 

For the massive objects, BUCnc and BNCnc, the predicted 
unstable ranges are compatible with the observed results (+ln), 
for both ttNL = 1-50 and = 1.89. The results for both objects 
are thus not sensitive to the value of otnl. This is to be expected 
because these more massive stars have shallower outer convection 
zones and thus their structures are less sensitive to the assumed 
value of the mixing-length parameter. In contrast, unstable ranges 
predicted for BW Cnc & BS Cnc show a strong dependence on this 
parameter. The discrimination between the two Cnl values reported 
in M99 is then partially reproduced. Finally, for the less massive 
object, B V Cnc, the values for q-nl cannot be distinguished. Never- 
theless, the observed ranges agree with the theoretical predictions 
within ±\n. Therefore, in general, these results constitute a consis- 
tent solution in terms of physics and cluster membership, and the 
observed and theoretical ranges of radial orders are in reasonable 
agreement for all the stars considered in this work. 

At this point, we should keep in mind that, as in M99, the 
ranges of unstable modes were estimated from non-rotating models 
for the observed stars. Therefore, we may wonder how accurately 
the so obtained theoretical ranges of unstable modes represent the 
observations of the rotating stars. However, the lack of theories de- 



scribing the effect of rotation on mode stability makes it difficult to 
overcome this problem. Nevertheless, assuming that mode stability 
depends predominantly on the effective temperature of the models 
(Pamiatnvkh 1975), it is possible to make a rough estimate of the 
elTect of varying the model's effective temperature on the location 
of the range of unstable radial modes. The effective temperature 
varies if the photospheric parameters are corrected for the effect of 
rotation (see Section [X2t . Keeping the remaining model parame- 
ters (mass, luminosity and chemical composition) constant we re- 
calculate for each model the linear growth rates. The so obtained 
ranges of unstable radial modes are compared with the observa- 
tions in Fig.O 

In general, the results are similar to those in Fig. [3J3, but the 
ranges of unstable modes are shifted. For a^i^ = 1.50 the unsta- 
ble ranges are shifted to lower values in the case of BU Cnc and 
BNCnc, but remain almost unaltered for the remaining stars. On 
the other hand, for = 1-89 the unstable ranges are shifted to 
higher values of n for the more massive objects, and to lower val- 
ues for the rest of stars. In particular, the discrimination pattern 
of BWCnc and BS Cnc in Fig.jSj) (obtained from the non-rotating 
equilibrium model parameters) is similar to the pattern of BU Cnc 
and BNCnc in Fig. |4] (obtained from model parameters for which 
Teff is corrected for rotational effects), which are more massive and 
hotter stars (see Table [3)- The effective temperatures of BW Cnc 
and BS Cnc in Fig. |4] are similar to the effective temperature of 
BVCnc in Fig.|3]and, as expected, the stability results are compa- 
rable. 



6 CONCLUSION 

In this work we studied ranges of unstable modes predicted by a 
stability analysis and compared the results with observations for a 
selected sample of 6 Scuti stars in the Praesepe cluster. With im- 
proved descriptions for rotational effects in both the evolutionary 
and adiabatic oscillations computations, we searched for a single 
consistent solution in terms of physics and cluster membership in 
order to match reasonably the observed and theoretically predicted 
ranges of overstable low-degree modes for all considered stars. 

Consistent solutions were found for models with an age of 
650 Myr, and which were computed with a mixing-length param- 
eter of Q-MLT = 1.614 and with a value for dov that is compatible 
with the grid of solutions provided by M99. In addition, for stel- 
lar models with log Teff = 3.87-3.88 a value of q-nl - 1-50 leads 
to a reasonable agreement between theoretical predictions and ob- 
servations of ranges of unstable modes, indicating that for these 
stars a smaller mixing-length parameter q-nl is required than sug- 
gested from a calibrated solar model. The need of a smaller value 
fo r g^L than that from a calibrated sola r model was also reported 
bv lPaszviiska-Daszkiewicz et al.l dlOOSl) for the d Scuti star FG Vir. 

Existing stability studies (including the present work) rely on 
calculations of non-adiabatic quantities without the inclusion of 
rotational effects; considering these effects in future investigations 
would lead to a substantial improvement of the stability analysis. In 
a first approximation we noticed an effect on mode stability through 
the change in effective temperature if it is corrected for rotational 
effects (see Section [T2t . The results presented here (Section [5. 2| > 
suggest, at least for some objects, a significant shift in the pre- 
dicted ranges of unstable modes. Conclusions going beyond these 
findings, however, would require the development of a theory that 
describes the interaction between rotation and pulsation dynamics. 
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